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Abstract 

Photo and DIS production of J/\E' are investigated and compared with calculations based 
on pQCD in the LLA approximation without and with screening corrections. The cal- 
culation includes corrections induced by the real part of the production amplitude, the 
skewed (off diagonal) gluon distribution function and the relativistic Fermi motion within 
the charmonium system. Our pQCD results are also compared with the predictions ob- 
tained from a Regge type two Pomeron model. Our results show that the screened pQCD 
model gives a better reproduction of the data than the non screened model. However, the 
available data does not enable us to exclude any of the three models we have examined. 
The predictions of these models, when extrapolated to both low and very high energies 
may provide a more effective discrimination between the different parameterizations. 
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1 Introduction 



The measurement of high energy photo and deep inelastic scattering (DIS) exclusive production 
of J/\& serves as an important testing ground for checking the validity of pQCD in the limit of 
a relatively small hardness and small x. The hardness of this process is fixed by m 2 ~ |M|/^, 
consequently |l|] HQ it is a suitable means of investigating the applicability and possible need 
for a re-formulation of hard pQCD when approaching the kinematic interface with the less 
understood npQCD soft domain. Indeed, over the past few years we have been witness to 
vigorous experimental, phenomenological and theoretical investigations of this process, which 
supplement the detailed analysis of F^x, Q 2 ), the proton structure function, and its logarithmic 
derivatives in the small Q 2 and small x limits. 

In the following we present a detailed study of J/ty photo and DIS exclusive production. 
Being relatively well measured, it may help to check the self consistency of, and discriminate 
between, relevant theoretical models . We recall that a key ingredient for determining the 
parton distribution functions (pdf) is the pQCD analysis of inclusive DIS which utilizes the 
DGLAP evolution equations^]. These pdfs are then used as input for the pQCD calculation of 
exclusive DIS channels (such as J/^> production), which are usually executed in the color dipole 
approximation^]. We shall mainly be interested in the correlated problems of determining the 
gluon saturation scale, the role of screening corrections (SC) and the relative importance of the 
soft npQCD component at the kinematic edge where the hard pQCD is applicable. 

As we shall specify in the next section, a realistic calculation of J/ty photo and DIS production 
depends on a few corrections to the original pQCD estimate: 

1) K\ accounting for the real part of the production amplitude @. The pQCD calculation of 
Eq. (|j) below relates to the imaginary part only. 

2) Kq denoting the correction resulting from the contributions of the off diagonal (skewed) 
gluon distributions 0. 

3) Kp denoting the correction due to Fermi motion of the heavy quark within the quarkonium 
system ||. As such, it accounts for the deviations from the simple static non relativistic estimate 
of the vector meson wave function. 



These corrections contribute significantly to the final (amplitude squared) estimate. Thus, in 
our investigation we shall assess possible ambiguities in our calculations which may be associated 
with the estimated corrections to the bare pQCD calculation. Our analysis follows the same 
lines as those of our recent paper || on <9F 2 /<9 In Q 2 • 



2 Photo and DIS production of J/ty 

The procedure for calculating the forward differential cross section for photo and DIS production 
of a heavy vector meson in the color dipole approximation is straightforward. The calculation is 
performed HH P] [P| ||10[ | II| in the LLA, assuming the produced vector meson quarkonium system 
to be non relativistic. The contribution of pQCD to the imaginary part of the t — differential 
cross section of photo and DIS production of heavy vector mesons is given by 



i 



where, xG DGLAP is the gluon distribution function as obtained from the DGLAP analysis. In 
the non relativistic limit we have 
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In the following we discuss the photo and DIS production of J/ty as there is an abundance 
of data available in this channel fll2f ||13|| [fP|] spanning a relatively wide energy range. From 



a theoretical point of view, its hardness is comparable to those we have investigated in our 
dFij 'd\\yQ 2 analysis || . In this study we have shown that the recent HERA data on the Q 2 
logarithmic slope of F 2 is well reproduced by DGLAP with either the CTEQ5HQ pdf input p~5j 



or by the GRV98NLO input [0 which is corrected for SC[17|. As we wish to maintain the 



compatibility of the F 2 and the J/\& production interpretations, we have confined our pQCD 
calculations only to the above two pdfs. In order to compare with the experimental data, 
which are given as integrated cross sections, to Eq. ([!]), we need to know B - the J/ty forward 
differential cross section slope. The experimental values are approximately constant with a 
possible moderate energy dependence [|12|]. Theoretically, each of the models we shall consider 
has a somewhat different estimate of B which we shall specify. 

The main difficulty with a pQCD analysis of J/^> is the observation that the simple dipole 
calculation needs to be corrected for the following reasons: 

1) A correction for the contribution of the real part of the production amplitude. This correction 
is well understood and is given by 

Kl = (1 +P 2 ), 
p = ReA/ImA = tg(*£), 

A = d\n(xG DGLAP )/d\n(±). (3) 
We note that A has a mild dependence on x at a fixed Q 2 . 

2) A correction for the contribution of the skewed (off diagonal) gluon distributions. This 
correction is calculated^] to be 

2 / 2^r(A + 2.5) \ 2 
Kg = { v^r(A + 4) ) ■ (4) 

3) A more controversial issue relates to the non relativistic approximation assumed for the 
J/ty charmonium. Relativistic effects, produced by the Fermi motion of the bound quarks, 
result in a considerable reduction of the calculated pQCD cross section||. The correction, 
Kp, is very sensitive to the value of m c . Ref. assumes that m c ~ 1.50 GeV and obtains 
K 2 F ~ 0.25 with almost no energy dependence. However, we note that the calculation of K F , 
regardless of its detailed construction, is very sensitive to the c-quark mass since a small change 
in the input value of m c changes the estimate of Kp significantly. Clearly, if we assume that 
m c = \Mj/q, ~ 1.55 GeV, i.e. a change of only about 50 MeV relative to the value assumed 
in Ref. ||, the Fermi correcting factor is identical to 1. We will, therefore, consider Kp as a 
free (energy independent) parameter. As we shall see, a calculation based on CTEQ5 gives 
Kp ~ 1.00, whereas a calculation based on GRV98 with SC gives Kp = 0.70. This value 
corresponds to a c-quark mass of approximately 1.52 GeV. 
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Based on the above, the expression for the integrated cross section is written 



2 2 2 1 fda\ pQCD 
a{i'p^J/9p) = K%.Kt l -K%---^j ^ , (5) 

/ . \pQCD 

where is given by Eq. (hj), B (the forward differential slope) is taken from the data 

and Kp is a free parameter. 



We wish to present an analysis for J/ty photo and DIS pQCD cross sections which is compatible 
with our recent investigation |J of the Q 2 logarithmic slope of F 2 . Note that the F 2 DGLAP 
analysis has been performed in NLO while the present calculation is carried out in the LLA 
modified by the above corrections. In the F 2 investigation we have shown [^] that the recent 
logarithmic slope data is well reproduced by three, rather different, formulations: 

1) A DGLAP calculation with CTEQ5 pdf input. This calculation has no explicit soft contri- 
bution. 



2) A DGLAP calculation with GRV98 pdf input modified by SC which are calculated jn] in 
the DLA. This calculation, as well, has no explicit soft contribution. 



3) A Regge type two Pomeron parameterization ||18|| in which the hard contribution is provided 



by a hard Pomeron, whereas the soft contribution is presented by a soft Pomeron. The two 
trajectories are parametrized to be 

a H {t) = 1.44+ 0.1* 

a s (t) = 1.08+ 0.25t (6) 



In the following we check and compare the above models with the J/\I/ cross sections. Our 
data base has 49 photo production points, out of which 35 points are HERA data. We have 
also studied just the 30 points at the high energy end with W > 50 GeV. The DIS data has 47 
data points, all of which come from HERA and have W > 50 GeV. 

1) Our non screened pQCD calculation, denoted CTEQ5NSC, is based on Eq. (|5p with CTEQ5 
input for xG DGLAP (x,Q 2 ). Our calculations for photo production, compared with the experi- 
mental data, are presented in Fig.l. This calculation has only a hard sector and as such in our 
calculations we took a fixed B = 4.73 CreV -2 [|I3J. We have adjusted the free parameter Kp so 
as to fit the data base of the integrated J/ty photo production and DIS cross section points. 
The best fit has no Fermi suppression, i.e. Kp = 1.00, and has the following values for the 
corresponding 

2 

a) For the complete data base ^ = 2.12. 

2 

b) For the photo production data ^ = 3.03. 

2 

c) For the HERA photo production data ^ = 1.14. 

2 

d) For high energy photo production (W > 50 GeV") ^ = 0.92. 

As seen in Fig.l, CTEQ5NSC overestimates the low energy data. In general, the overall x 
dependence of CTEQNSC is somewhat softer than the harder x behaviour suggested by the 
data. We conclude that CTEQ5NSC does not produce a good fit to the complete data base. 

2 

However, its ^- is considerably improved once we ignore the low energy data. 



3 



2) For the SC calculation, done in the DLA, we follow our earlier publications |19| and define 



the damping factors due to the screening in the quark sector i.e. the percolation of the cc 
through the target. This is given by the following expressions for the longitudinal and trans- 
verse damping factors 



= | — ^— I (?) 



and 



1 + (i - D E ^y 



D 2 qT =\ ^ * |, ($) 



where, for N c = 3, we have 



^ = ^m xGDGLAP ^ Q2) - (9) 

The above expressions are derived[|ll| assuming that 7 << 1, where 7 is the DGLAP anoma- 
lous dimension. As a result of this approximation we overestimate the SC in the exceedingly 
small x limit. However, within the kinematical range of this investigation this excess is small 
enough to be neglected for x > 10~ 3 . The correction for smaller x may be as large as 10%. 
Actually, the exact calculation makes our results moderately harder and consequently improves 
our reproduction of the data (see Fig. 2). We do not consider this to be important as it will 

2 

result in a small reduction of our present ^ values, which are excellent as they stand. Our 
expression for D 2 , the damping in the gluon sector, is the square of the gluon damping used in 
Refs.P][p|. Our final expression for the integrated cross section is 



pQCD 

a(j*p^J/*p) = K\-Kl ±— .(^ ] -Dl-Dt (10) 



'day 
B(R 2 ) \dt) t __ 



where D q denotes the L and T components as appropriate. As in our F2 study Hfl?]], for 
xG DGLAP (x, Q 2 ) we use the GRV98 pdf input. 

Our calculations as compared with the data are presented in Figs. 2 and 3. with an adjusted 
value of K F = 0.70. We get the following values for the corresponding 

2 

a) For the complete data base ^ = 0.94. 

2 

b) For the photo production data ^ = 0.91. 

c) For the HERA photo production data ^ = 0.87. 

2 

d) For high energy photo production (W > 50 GeV) ^ = 0.58. 

As is visible from Fig.2, the suppression induced by the SC is appreciable even though n g < 1, 
i.e. below gluon saturation. This is consistent with the general observation ||20|| that SC, which 



are the consequence of high gluon density, precede the gluon saturation state. Note that 
R 2 = 8.5GeV~ 2 , which is the essential parameter in the SC calculation presented both in 
Ref. M and here, is not a free parameter but is determined directly from the J/^ photo pro- 
duction forward slope. In a SC model, such as ours, we expect J2~H| a weak dependence of B on 
x. This is demonstrated in Fig. 4 together with the relevant HERA data. For a non screened 
calculation B is fixed at \R 2 - 
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J/^ Photoproduction 





Figure 1: Photo production of J/^f as a function of W andx. Data and CTEQ5NSC calcula- 
tions. 
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J/^ Photoproduction 




Figure 2: Photo production of J/ty as a function of W andx. Data and GRV98 calculations 
with and without SC. 
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Figure 3: DIS production of J/^. Data and our calculations. 
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3) Donnachie and Landshoff (DL) have followed their F 2 analysis|18| with an updated pub- 
lication on charm production |22| in which they show that the J/ty photo and DIS data on 
the integrated and differential cross sections are compatible with their model, which we de- 
note DL2P. We discuss the differences between DL2P and the predictions of CTEQ5NSC and 
GRV98SC in the next section. 



3 Discussion and Conclusions 

The results of the present investigation corroborate the conclusions obtained in our previous 
study of dF 2 jd In Q 2 , where we have shown that the available data is consistent with either 
CTEQ5NSC, GRV98SC or DL2P in the kinematic range of Q 2 > 1.9 GeV 2 and x < 1CT 2 . A 
comparison of the photo production (Q 2 = 2AGeV 2 ) predictions of these models over a wide 
x range is presented in Fig. 5. Following are some concluding remarks 

1) The differences between CTEQ5NSC and GRV98SC, within the available experimental kine- 
matic window, are relatively small. We note, though, a systematic difference at low energies 
{x > 10~ 2 ) where CTEQ5NSC is considerably larger than GRV98SC. This difference is respon- 
sible for the high ^ values we found for CTEQ5NSC. Similar differences also exist for the Q 2 

logarithmic slope predictions of these models || in the low energy (high x) limit, but relevant 
data is not available. From the study of J/^f photo production we conclude that the GRV98SC 
parameterization of this process extrapolates well to the very low energy domain in which the 
soft (npQCD) sector plays an increasingly important role. On the other hand, CTEQ5NSC 
does not have this property and it overestimates the low energy data. 

2) We also note a difference between CTEQ5NSC and GRV98SC in the high energy region 
(x < 10~ 4 ). This was not observed in our predictions for the logarithmic F 2 slope and is at- 
tributed to our approximate calculation of D 2 (J/ty) which results in a small excess of the SC 
for very small x. As we noted earlier, an exact calculation of D 2 (J/^) will improve the results 
of GRV98SC and bring it close to those of CTEQ5NSC. 

3) As we have shown, CTEQ5NSC provides an acceptable reproduction of the high energy 
photo production data. These results are obtained with K F = 1.00 which implies that the J r /\& 
charmonium bound system is strictly non relativistic, or else, that the normalization of the 
calculation requires some adjustment. We have further investigated this model by applying SC 
to the calculation with Kp = 1.00, denoting it CTEQ5SC. Our output is non satisfactory as 
it systematically underestimates the data. This result supports our suggestion that CTEQ5 
may contain significant screening effects which are absent in the boundary conditions used in 
GRV98. In this context, it is suggestive to assume that CTEQ5SC results are too small due to a 
possible double counting of the SC. Our calculated SC introduce a 20% deficiency at x = 10~ 4 . 
We consider this as the uncertainty in the CTEQ5 parameterization in the exceedingly small x 
limit. This is compatible with the similar conclusion reached in Ref.@. The pQCD calculation 
used in this investigation is performed in the LLA and is, thus, subject to corrections. We 
distinguish between the normalization of the models we have discussed, which may change, and 
their hardness as reflected in the energy, or x, dependence, which is a more stable property. 

4) Within the available kinematic photo production window, DL2P predictions which differ 
somewhat from the results of the two pQCD models, still provide an adequate reproduction 



of the data||22|| . At this stage, with the given experimental error bars, one cannot conclu- 
sively exclude any of these models. At higher energies DL2P is significantly larger than either 
CTEQ5NSC or GRV98SC. This feature of DL2P is a consequence of the hard Pomeron com- 
ponent which has a very high Regge trajectory intercept. As a result its energy dependence is 
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Figure 4: The energy dependence of the forward differential slope of J/ty photo production. 
Data and SC calculations with several values of R 2 . 
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much harder than the pQCD models whose A values are considerably smaller. This is clearly 
seen at x — 10 -5 where the DL prediction is about 5 times larger than that of GRV98SC. This 
x range will become accessible at THERA and, thus, enable an experimental discrimination 
between DL2P and the pQCD models. The same behaviour of DL2P has also been observed|9|] 
in the small x limit of di^/dlnQ 2 . 

5) We have investigated the possible role of a soft component, parameterized in the DL form, 
which was added to GRV98SC. Our analysis shows that such an addition does not improve the 

2 

We, thus, conclude that any soft component in our model is exceedingly small. 

6) The three models discussed have different predictions for B, the J/ty forward differential 
slope. CTEQ5NSC is a model with only a hard, non modified, sector. As such, we expect its B 
to be a constant for which we took the experimental value. GRV98SC is also a model with only 
a hard sector but it is corrected by SC. As a result, the model predicts a modest shrinkage of B 
which is presented together with the data in Fig.4. DL2P is a sum of a soft component, which 
dominates in the low energy limit and a hard component which dominates in the high energy 
limit. Since the two Pomerons have different trajectory slopes, DL2P predicts an anti shrinkage 
effect, where B decreases monotonically with W. The recent ZEUS|12| results contradict this 
feature. 



Overall, our investigation shows that GRV98SC results in a better reproduction of the data 
than the other two models considered. However, the available data does not enable us, as yet, 
to conclusively exclude any of the models we have examined. 
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